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We report results from Monte Carlo simulations of the two- and three-dimensional gauge glass at 
low temperature using parallel tempering Monte Carlo. In two dimensions, we find strong evidence 
for a zero-temperature transition. By means of finite-size scaling, we determine the stiffness exponent 
6 = —0.39 ± 0.03. In three dimensions, where a finite-temperature transition is well established, we 
find — 0.27±0.0I, compatible with recent results from domain-wall renormalization group studies. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

The gauge glass is a model often used to describe the 
vortex glass transition in high-temperature superconduc- 
tors. Still, there are areas which need to be understood. 
In two dimensions, there is an ongoing controversy as to 
whether a spin-glass transition occurs at finite temper- 
ature or not. In three dimensions, a finite-temperature 
transition is well established^, although there is no con- 
sensus on the value of the stiffness exponent. 

Recently, claims of evidence for a finite-temperature 
transition in two dimensions with Tc = 0.22 by Kim^ us- 
ing resistively shunted junction (RSJ) dynamics, and by 
Choi and Park^ who study the scaling of the spin-glass 
susceptibility via Monte Carlo simulations, have been 
made. In contrast, Granato^ and Hyman et al^, who 
also use RSJ dynamics find evidence of a zero-T transi- 
tion. In addition, Monte Carlo simulations by Fisher et 
al.^ and Reger and Young^ show evidence of a T = 
transition, although the simulations were not performed 
at low enough temperatures. 

Here we perform Monte Carlo simulations of the 
two- and three-dimensional gauge glass, using parallel 
temperingSiS, to go to significantly lower temperatures 
than was possible in earlier work. In particular, we 
are now able to cover the temperature range in two di- 
mensions where the claimed spin-glass transitionSi^ takes 
place. We find strong evidence that Tc = 0. 

In three dimensions, we study the gauge glass at very 
low but finite temperatures to provide a good estimate 
of the stiffness exponent 9. Earlier estimates, which were 
obtained from ground-state methods, are inconsistent. 
One groupiSiiiiiSii^ finds values consistent with « 0, 
whereas anotherii^ii^*iSiii finds 9 in the range 0.26 - 
0.31. We find that 9 = 0.27±0_01 which agrees with the 
results of Refs.H and ll4ll5ll(il7i 

II. MODEL, OBSERVABLES, AND METHOD 

The Hamiltonian of the gauge glass is given by 

H ^ -^cos{(j)i - - Aij), (1) 



where the sum ranges over nearest neighbors on a hyper- 
cubic lattice in D dimensions of size A'^ = , and (f>i 
represent the angles of the XY spins. Periodic boundary 
conditions are applied. The Aij are quenched random 
variables uniformly distributed between [0, 27r}i^. Be- 
cause Aij represent the line integral of the vector po- 
tential between sites i and j, we have the constraint that 
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Traditionally, one uses the Binder ratic^^ to estimate 
the critical temperature Tc- As for the gauge glass the 
Binder ratio cannot exceed unitjii, the splaying of the 
curves is small, and Tc is difficult to establish. In order 
to avoid this problem, we use a method introduced by 
Reger and Youngifi in which one calculates the current / 
defined as the derivative of the free energy with respect 
to an infinitesimal twist to the boundaries, i.e., 

I{L) = —'^sm{(t)i - (jji+x ~ Aii+x) . (2) 

i 

In this case, the twist is applied along the x-direction. 
As [(/(L))]av — 0, we actually calculate the root-mean- 
square current /mis- Here (■ • ■) represents a thermal av- 
erage, whereas [• ■ -Jav represents a disorder average. Note 
that the current scales^i as 

Irn.s=iiL'/''{T-Tc)) , (3) 

for Tc> 0, whereas for Tc = 0, 

= L-^'''i{L^I''T) . (4) 

Equation Q indicates that the T = stiffness exponent 
9 is negative and equal to —1/v (as /mis ^ L^)- Equation 
Q shows that if Tc > 0, the curves for /j-ms for different 
sizes intersect at the critical point, whereas Eq. I^J shows 
that if Tc = the data decrease with increasing size at 
T = 0. For a finite-temperature transition we expect 
6* > 0, since then the ordered state at T = is "stiff" on 
large scales, and so will presumably resist small thermal 
fluctuations. On the other hand, for a zero-temperature 
transition, the system will then easily break up under the 
influence of thermal fluctuations and therefore 6* < 0. 

We also have calculated the spin-glass susceptibility, 
XsG ' defined by 

xso = mq^)U , (5) 
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where q is the spin-glass order parameter: 



1 ^ 
N ^ 



exp[z(0f - 4)\ . 



(6) 



Here, a, (3 are two rephcas of the system with the same 
disorder. According to standard finite-size scahng the 
spin-glass susceptibility, defined in Eq. 0, behaves as 



(7) 



meaning that at criticality (T — Tc) it diverges with a 
power law, i.e., Xsg ^ L'^^^ . Note that the power-law 
prefactor in Eq. Q with an unknown exponent compli- 
cates the analysis of Xsg compared with that for /rms- 

For the simulations, we use parallel tempering Monte 
Carlc&fi. as it allows us to study systems at lower tem- 
peratures than with conventional methods. For details 
about the implementation, equilibration tests, and de- 
tailed parameters of the simulations we refer the reader 
to Ref. 20. In two dimensions, the highest temperature is 
1.058, whereas the lowest temperature is 0.13 for L = 4, 
6, 8, 12, and 16 (for L = 24 the lowest temperature is 
0.20). The number of temperatures Nt = 30 is chosen 
to give satisfactory acceptance ratios for the Monte Carlo 
moves between the temperatures {Nt = 24 for L = 24). 
In 3D the lowest temperature studied is 0.05 (note that 
Tc « 0.45)i whereas the highest temperature is 0.947, 
and Nt = 53 for L = 3, 4, 5, 6, and 8. 



III. RESULTS FOR D = 2 

Figure ^ shows a finite-size scaling plot of the data 
for /rms in two dimensions according to Eq. 10}. The 
plot shows that the data collapse if L^^^T is small for 
all values of L, and over the whole range of L^I^'T for 
the largest sizes. The inset shows results of the unsealed 
data. At all temperatures, the data decrease with in- 
creasing L indicating, from Eq. (O, that Tc must be less 
than the range of temperatures studied. The data in 
Fig. n are therefore consistent with a zero-temperature 
transition, but with significant corrections to scaling at 
intermediate temperatures. We estimate the stiffness ex- 
ponent to be 6* = —0.39 ± 0.03. The above error bar is 
estimated by varying slightly until the data do not col- 
lapse well. This result is consistent with recent work of 
Akino and KosterlitaM who find 9 = -0.36 ± 0.01. Choi 
and Park's^ parameters [T^ = 0.22, l/v = 0.88) yield 
very poor scaling, especially near the proposed Tc. In 
fact, we are unable to get a reasonable fit to the data for 
/r,ns according to Eq. Q for any finite Tc- 

Figure 121 shows a scaling plot of the data for Xsg 
cording to Eq. iQ for Tc = and l/v = 0.39, the same 
parameters found in the scaling of /rms, together with 
T] = 0.0 ± 0.1, which is expected at a zero-temperature 
transition. The data at low T and for the largest sizes 
scale well, but the data away from this range show devi- 
ations. Allowing 1/v to vary, we find 1/v = 0.50 ± 0.03. 
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FIG. 1: Scaling of the root-mean-square current /mis in two 
dimensions according to the form expected if Tc = 0, Eq. 
We see acceptable scaling of the data at low temperatures. 
Deviations at higher T are presumably due to corrections 
to scaling. This plot is for 9 = —Ijv — —0.39. The inset 
shows 7rms as a function of T for different system sizes. At all 
temperatures, the data decrease with increasing L indicating, 
from Eq. JHJ, that if Tc is finite it must be less than 0.13. 
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FIG. 2: Scaling of Xsg according to Eq. Q with Tc = 0. 
The data for large sizes and low temperatures collapse with 
77 = and l/v — 0.39. The inset shows a scaling plot for the 
optimal value 1/1/ = 0.50 (and -q = 0). For these values of 
exponents, the collapse extends to a larger range of sizes. 



The inset shows data for this optimal value, where only 
the L — 4 and 6 data are not part of the scaling function 
for all L^^"T. The fact that the best values oi 1/v are 
not precisely the same when obtained from Xsg ^'^d /rms 
presumably indicates that scaling is only valid for fairly 
low temperatures and large sizes, and that, despite our 
working at quite low temperatures, we have only a lim- 
ited range of data which are fully in the scaling regime. 
We also scale the data for Xkg with Choi and Park's 
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FIG. 3: Effective stiffness exponent 9cs{T) as a function of 
temperature in three dimensions for low temperatures. The 
data extrapolate linearly to T = 0. 
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parameters^. The data collapse is poor near Tc- The 
scaling of /rms is also niuch worse, indicating that /rms 
distinguishes between a finite Tc and Tc = much better 
than XsG because of its simpler finite-size scaling form. 
An attempt to scale the data for Xsg with Tc = 0.13, the 
lowest temperature simulated, shows that the best fit is 
obtained with = 0.68 and rj = 0.19. While the data 
scale acceptably well, the data for Irms scale poorly. 

According to Eq. Q for T — Tc, the data for Xsa 
should lie on a straight line at Tc- However, the data 
in the vicinity of T = 0.22, the transition temperature 
claimed by Kim^ and Choi and Park,'' show a strong 
downward curvatureSS in our analysis, indicating that 
this is actually above Tc- Only around the lowest tem- 
perature where we have data, T = 0.13, is the curvature 
small, although it still greatly exceeds the error bars. 



This indicates that Tc < 0.13, which is compatible with 
our data for Irms- 



IV. RESULTS FOR D = 3 

Olson and Young^ have investigated the critical region 
of the three-dimensional gauge glass obtaining a lower 
bound for the stiffness exponent oi 9 > 0.18. Some pre- 
vious resultsi2iiiii2J^ find in the range < 9 < 0.077 
whereas othersi^iiSiiSiil find a much larger value, 0.26 < 
9 < 0.31. 

We can estimate 9 from our data for Ij-ms since I^ms ^ 
when L^/'^{T — Tc), the argument of the scaling func- 
tion in Eq. Q, tends toward infinity. To obtain an esti- 
mate of 9, we perform a linear least-squares fit of ln(/i.ms) 
against ln(L) for each temperature in order to obtain an 
effective stiffness exponent 9es{T) which depends on the 
temperature. Figure |21 shows that 9ca{T) can be fitted 
well to a linear form at low temperatures. Extrapolating 
to T = 0, we obtain 9 = 0.27±0.01, which is clearl y p osi- 
tive and consistent with the results of Refs.lll ll4iri5Bl6lll7l 



V. CONCLUSIONS 

Our results from Monte Carlo simulations of the two- 
dimensional gauge glass are consistent with a T = 
transition with a stiffness exponent 9 = —0.39 ± 0.03. 
These results are incompatible with the prediction made 
by Kim.2 and Choi and Parl^ that Tc = 0.22. In three 
dimensions we report the first reliable estimate of the 
stiffness exponent from finite-temperature Monte Carlo 
simulations. We find 9 = 0.27 ± 0.01, which agrees with 
the results of Refs. 
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